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ABSTRACT 

With observations from the Interface Region Imaging Spectrograph (IRIS), we track the com¬ 
plete evolution of ~11 MK evaporation flows in an Ml.l flare on 2014 September 6 and an XI. 6 
flare on 2014 September 10. These hot flows, as indicated by the blueshifted Fe xxi 1354.08A line, 
evolve smoothly with a velocity decreasing exponentially from ^200 km s -1 to almost station¬ 
ary within a few minutes. We find a good correlation between the flow velocity and energy 
deposition rate as represented by the hard X-Ray flux observed with the Reuven Ramaty High 
Energy Solar Spectroscopic Imager (RHESSI), or time derivative of the soft X-Ray flux observed 
with the Geostationary Operational Environmental Satellites (GOES) and the HINODE X-ray 
Telescope (XRT), which is in general agreement with models of nonthermal electron heating. 
The maximum blue shift of Fe xxi appears approximately at the same time as or slightly after 
the impulsive enhancement of the ultraviolet continuum and the Mg II 2798.8 A line emission, 
demonstrating that the evaporation flow is closely related to heating of the lower chromosphere. 
Finally, while the hot Fe xxi 1354.08 A line is entirely blueshifted with no obvious rest compo¬ 
nent, cool chromospheric and transition region lines like Si IV 1402.77A are often not entirely 
redshifted but just reveal an obvious red wing enhancement at the ribbons, suggesting that the 
speed of chromospheric condensation might be larger than previously thought. 

Subject headings: Sun: flares—Sun: chromosphere—Sun: transition region—line: profiles—magnetic 
reconnection 


Introduction 


Solar flares are among the most intense energy 
release events in t he solar atmosphere. In the stan¬ 
dard flare mode l (I Car michaell 1964: Stu rrockll19661: 
Hiravam5ll974 IKopp fc Pneuman Il976 h energy 


released from magnetic reconnection first heats 
the plasma and accelerates particles in the corona. 
The released energy is then transported along the 
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newly reconnected flare loops towards the chro¬ 
mosphere through either thermal conduction or 
nonthermal electrons, which results in an intense 
and rapid heating of the chromospheric material 
up to a temperature on the order of 10 MK. The 
chromospheric heating is revealed through bright 
flare ribbons in images of Ha and UV contin¬ 
uum. At the same time the nonthermal electrons 
can produce hard X-Ray (HXR) emission through 
bremsstrahlung with the thermal ions. The over¬ 
pressure associated with the chromospheric heat¬ 
ing drives the heated plasma upward to fill the 
flare loops. This so-called ” chromospheric evapo¬ 
ration process” has been proposed to explain the 
soft X-Ray (SXR) emission which usually peaks 
after the HXR peak (e.g., lNeupertlll968 ). In the 
impulsive phase of some flares the time history 
of HXR flux has been found to closely match the 
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derivative of SXR light curve. This relationship is 
called the Neupert effect and is often quoted as a 
manifes tatio n of ch romospheric! evaporation (e.g., 


manuestation ot cn romospnenc evaporatio n (e.g. 

Hudsonl 1991 : Dennis fe Zarrolll993l : lVeronig et al 


2005; Liu et al. 2' 


There have been only a few broad-band imag- 
ing observations of chromospheric evaporation 


(Silva et al 

1997 

:lLiu et al. 2006; Ning et al. 20091 

Nitta et al. 

2012 

; Zhang & Ji 20131). These au- 


thors found X-Ray sources moving towards the 
loop top during the impulsive phase of the flares 
and interpreted them as thermal emission at the 
evaporation front or nonthermal emission due 
to decreasing stopping distance from the loop 
top caused by the increasing density. Most of 
our knowledge about chromospheric evaporation 
is based on spectroscopic observations, proba¬ 
bly because spectroscopy is the best way to ob¬ 
serve vertical flows in on-disk observations. The 
most direct spectral signature of the evapora¬ 
tion how should be the blue shift of hot emis¬ 
sion lines. Indeed, early observations by the X- 
Ray spectrometers mounted on the P78-1, the 
Solar Maximum Mission and the Yohkoh space¬ 
craft have detected blueshifted components of 


25 MK (|Doschek et al. 


; Feldman et al. 1980t 

Antonucci et al. 19821: 

Mariska & Doschek 19931). 


indicating the presence of hot and fast (up to 
^400 km s -1 ) plasma uphows. Enhancement in 
the blue wing of the Fe xxi 1354.08A line in¬ 
dicating an uphowing component at a speed of 
~200 km s -1 was also detected by the Ultraviolet 
Spectrometer and Polarimete r (UVS P) on board 
the Solar Maximum Mission ( Mason et al. 1986). 


In ah these observations, there is a dominant sta¬ 
tionary emission component which likely results 
from the surrounding stationary hot plasma in 
the hare loops. 

Observations from the Coronal Diagn ostic 
Spectrometer (CDS, lHarrison et al.l ll99 5) on 
board the Solar and Heliospheric Observatory 
(SOHO) have greatly improved our understand¬ 
ing of chromospheric evaporation. For example, 
iBrosiusI ( 200 3) detected blueshifted line prohles 


from the O ill, O V, Mg x and Fe xix ions early 
during the impulsive rise of a hare. He also found 
that ah lines except the Fe xix line subsequently 
exhibited multiple components, indicating down- 
hows resulting from cooling of the heated plasma. 


Harr a et al.l ( 2005 ) reported two ” ribbons” where 
the Fe xvi and Fe xix line prohles exhibited blue 
shifts of a few tens km s -1 . These two ” ribbons” 
were found to separate with time, which was in¬ 
terpreted as a signature of chromospheric evap¬ 
oration in sequentially reconnected loops. Blue 
shifts of some transition region lines have also 
been found before the impulsive phase of som e 
hares ( Brosius fe Phihip3l2Q04t lHarra et al.ll2QQ5h . 
Th e earl y evap o ratio n how has been proposed 
to suppo rt the breakout 


by lHarra et al 


mode l of coronal mass eruption ([Antiochos et al.l 
19991 ). where magneti c reconnection s tarts b efore 


the hlament eruption. iMilligan et al.l ( 2006a ) and 
Milligan et al.l ( 2006b[ ) presented exa mples of gen¬ 


tle an d explosive evaporation (e.g., iFisher et al 


Il985bl) in the impulsive phase of hares, respec¬ 
tively. In the former case, heating by relatively 
low hux of nonthermal electrons leads to evapo¬ 
rated plasma uphows at several tens km s -1 with 
no associated downhows. In the latter case, heat¬ 
ing by relatively high hux of electrons leads to 
an overpressure which results in both an upward 
expansion of the hot plasma at several hundred 
km s -1 (^230 km s -1 blue shift of Fe xix 592.23A) 
and a downward condensation of the cooler mate¬ 
rial at several tens km s -1 (^40 km_s_ -1 re d shift o f 
He I 584.33A and O V 629.73A). iDel Zanna et al" 


( 20061 ) also found that the Fe xix 592.23A line 
was largely blueshifted by ^140 km s -1 in an ex¬ 
plosive evaporation. The Fe xix blue shift was 


found to be r^j 80 km s 


-1 


in two h ares studied 
by Falchi et al.l (|2006f ) and Rafterv et al.l (|2009f ). 
The Fe xix 592.23 A line appears to be entirely 
blueshifted in many of these CDS observation. 


However, as pointed by Young et al. | (l2013h . this 


is likely because the low spectral resolution of CDS 
blurs the multiple components of the line. Indeed, 
Fe xix 592.23A line prohles with both a nearly sta¬ 
tionary component and a highly blueshifted (by 
~200 k m s -1 ) component were id e ntified in a few 
hares ( Teriaca et al.l 120031 120061: IMilligan et al.l 
2006b| ) . Coronal uphows at a speed of several tens 


of km s 1 have also been found by CDS_in the 


late gradual phase of a hare ([Czaykowska et al 


1999h . These uphows have been suggested to re¬ 


sult from chromosph eric evaporation driven by 
thermal conduction ( Czavkowska et al ]Hm|). 

With observations of t he EUV imaging spec¬ 
trometer (EIS, ICulhane et al.l 12007 ). our knowl- 
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edge about chromospheric evaporation has been 
further enhanced in the past 6 years. EIS has a 
better coverage of the coronal temperatures over 
CDS and allows study of the temperature depen¬ 
dent behavior of the evaporation flows. Although 
different patterns of evaporation flows have been 
found in different flares and even at different pix¬ 
els within the same ribbon, two surprising results 
appear to be reported for many flares. First, pro¬ 
files of hot lines from Fe xxm (formed at ~14 
MK) and Fe xxiv (formed at ~18 MK) usually 
exhibit a blueshifted component or blue wing en¬ 
hancement besides a nearly stationary compo¬ 
nent, (e.g., iMilligan fe Dennisl 120091 : lYoung et al. 
2013[ ), which is no t consistent with single fl are loop 


model s (e.g., Emslie fe Alexand er 1987; Lin et al 


2009at ) and may be explained by an ensemble of se- 
qnential heating of many s m all-scale threa d s (e.g 


Warren fe Doschekl 120051 : iReeves et al.l 12007 ). 


This result suggests that the spatial resolution of 
EIS (^3") is insufficient to separate the evapora¬ 
tion flow from the ambient stationary hot plasma. 
The velocities associated with the blueshifted 
components can be as high as ~400 km s _1 . 
In some flares blueshifted components are even 
present in profiles of warm lines such as Fe xm, 
Fe xiv, Fe xv and Fe xvi (formed at 1.5-3 


MK. Wat,a 

aabe et al. 2010: Del Zanr 

La et al. 2011; 

Li & Ding 

2011; Graham et al. 2011 

). The veloc- 


ities inferred from the profiles of these warm lines 
are mostly less than 150 km s -1 . In most cases 
the stationary component is dominant, although 
a dominant blueshifted co mponent has also been 
reported in a few flares ( Watana be et al.l 120101: 


Li fe Din 5 l201ll I Young et al.l 2013 ). We note, 


however, completely blueshifted Fe xxm line pro- 
files were identified occasionally in a few flares 


( Brosiusl 120131: iDoschek et al. 2013 ). Second, in 


explosive evaporations the temperature at which 
the Doppler shift turns from red to blue appears 
to be much higher than those predic ted b y some 
models. For example, the model of iFisher et al 


( 1985at ) predicts a condensation of plasma in the 
chromosphere. However, EIS observations often 
show a few tens km s -1 redshift of emission lines 
formed at m uch higher temperature s. For exam - 


Milligan fe Dennis 


and Young et al 


P le > T 

(120131 ) found that the transition occurs in the 
temperature range of 1.5-2.0 MK and 1.1-1.5 MK, 
respectively. A high transition temperature of 


^0.8 MK was repo rted even in a small B1.4 flare 
( Chen fe Ding|l2010j ). A recent model of Liu et al 


( 2009a ) shows that continuous energy deposition 
throughout the impulsive phase can heat the un¬ 
derlying chromosphere up to 2 MK, which may 
explain the observed red shifts of some coronal 
lines. Besides the Doppler shift, the nonthermal 
line broadening associated with the evaporation 
flow has also been investigated using EIS data. 
Significant nonthermal broadening (~100 km s -1 , 
FWHM) has b een found for the hot Fe xxm an d 
Fe xxiv l ines (Milliganllioill lYoung et al. 2013 ). 
iMilliganl ( 2011 ) also found a positive correlation 
between the blue shift and nonthermal width for 
the Fe xv and Fe xvi lines at one loop footpoint. 

With the launch of the I nterface Regio n Imag¬ 
ing Spectrograph (IRIS. iDe Pontieu et al.l 12014 ). 
chromospheric evaporation can now be observed 
at unprecedented high spatial (~0.33"), spectral 
(~3 km s -1 ) and temporal (^3 s) resolution. Al¬ 
though there is only one line, Fe xxi 1354.08A 
(for mation t emperature ~11 MK from CHIANTI, 
Landi et aD ( 2Q13h ). formed above 2 MK in the 


IRIS spectral range, significant new insight about 
chromospheric evaporation has already been re¬ 
vealed from IRIS observations. The most in¬ 
triguing finding is the completely blueshifted 
Fe xxi 1354.08A line profiles in the evapora¬ 
tion flows, which suggests that the resolution of 
IRIS is high enough to fully resolve the evap¬ 
oration flows. The reported blue shift of the 
Fe xxi 1354.08A line a ppears to be mostly around 
100 km s -1 or smaller ( Young et al. 2015; Li et al 


[2015; IPolito et, al.lEoT 
though a clear 
ported in a Cl.6 flare (Tian et al 


Sadvkov et al.l 12015 ). al 


blue_shift was re- 
2014al ). These 


subarcsecond-resolution observations often reveal 
a spatial offset between the continuum enhance¬ 
ment and the blueshifted Fe xxi emission, which 
is likely due to the fact that the hot evapora¬ 
tion fl ows are located higher in the loop legs. 


Polito et al. ( 2015 ) also found a decrease of both 


the blue shift and nonthermal line width of Fe xxi, 
although only six data points were measured due 
to the 75 s cadence of the raster scans. 

Here we report results from IRIS sit-and-stare 
observations of chromospheric evaporation in two 
flares. With a cadence of ^9.5 s, we track the 
complete evolution of hot (^11 MK) evapora¬ 
tion flows in the Ml.l flare on 2014 Septem- 
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ber 6 and the XI.6 flare on 2014 September 10. 
With simultaneous observations from the Reuven 
Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI. iLin et al. 200 2). the Geostationary Op¬ 
erational Environmental S atellites (GOES) and 
the X-ray Telescope fXRT. [Golub et aJ]l2007l ) on 
board Hinode, we also investigate the possible cor¬ 
relation of the flow velocity with the energy depo¬ 
sition rate and chromospheric heating. Our obser¬ 
vations reveal new insight into the chromospheric 
heating and evaporation processes during flares. 


2. The Ml.l flare on 2014 September 6 
2.1. Observations and data reduction 

An Ml.l flare occurred in NOAA AR12157 at 
about 17:00 UT on 2014 September 6. Two peaks 
are present in the GOES SXR light curves, sug¬ 
gesting two episodes of energy release in this flare. 
Here we focus on the first episode from roughly 
16:51 UT to 17:00 UT, when the plasma heating 
process was clearly observed by the IRIS slit at 
one footpoint of the flare loops. In this observa¬ 
tion IRIS pointed at (-732", -299"), with a roll 
angle of 45 degrees. The slit was fixed at the same 
location on the Sun throughout the whole observ¬ 
ing period (sit-and-stare, 11:23 UT - 19:30 UT). 
The spatial pixel size was 0.166". The spectral 
dispersion was ~0.026 A per pixel in the far ultra¬ 
violet (FUV) wavelength bands (FUV spectrally 
summed by 2) and ^0.025 A per pixel in the near 
ultraviolet (NUV) wavelength band. The cadence 
of the spectral observation was ^9.5 seconds, with 
an exposure time of 8 seconds. Slit-jaw images 
(SJI) were taken at a cadence of ^19 seconds in 
each of the 1400A and 1330A passbands, and the 
exposure time decreased from 8 seconds to 2.4 sec¬ 
onds after 16:55:15 UT. Dark current subtraction, 
flat field correction, geometrical and orbital vari¬ 
ation corrections have been applied in the level 2 
data used here. Several relatively strong neutral 
lines have been used for absolute wavelength cali- 
bration, e.g., O I 1355.59771, S I 1401.51361 and 
Mn i 2801.9081. These lines were assumed to be at 
rest outside the flare regions. The fiducial marks 
were used to achieve a coalignment of the two SJI 
passbands and different spectral windows. Fig¬ 
ure QJ A)-(D) shows the IRIS SJI 1330A image and 
spectral images in three wavelength windows at 
16:55:35 UT. 


The 131A passba nd of the Atmosph eric Imag¬ 
ing Assembly (AIA, iLemen et al.l l2012h on board 
the Solar Dynamics Observatory (SDO. BPesnell et al. 
2012h sam ples emission mainly from Fe XXI in 
flares (e.g., lO'Dwver et al.l i 2010). We use images 
of this passband for the examination of the flare 
loop morphology (Figure DJE)). We notice that 
early during the impulsive phase t his passban d is 



AIA images by 45 degrees to match the IRIS 
spacecraft orientation. We coaligned the AIA 
1600A (mainly FUV continuum and C iv) and 
IRIS 1330A (mainly FUV continuum and C il) 
images. AIA 131 A images were then automat¬ 
ically aligned with the IRIS images since AIA 
images in different passbands are automatically 
coaligned after applying the standard SolarSoft 
(SSW) routine aia_prep.pro. The cadence of the 
AIA 13lA images is 12 seconds. 

RHESSI data is also available for this flare. 
The RHESSI 25-50 keV image reconstructed with 
the CLEAN algorithm using detectors 2-8 from 
16:55:00 UT to 16:55:32 UT shows two HXR 
sources, which are clearly located at the two foot- 
points of the flare loops (Figure DjE)). Note that 
the RHESSI image has also been rotated by 45 
degrees. 

We present the time sequence of the IRIS ob¬ 
servation (including SJI 1330 A images and three 
spectral windows) and the corresponding AIA 
131 A images in an online movie (ml.mp4). The 
GOES SXR and RHESSI HXR light curves are 
shown in the movie. The RHESSI 25-50 keV 
HXR and the time derivative of GOES 1-8A SXR 
show a reasonably good match, indicating the 
existence of the Neupert effect. We also de¬ 
fine the NUV and FUV continua by averaging 
the IRIS spectra over slit positions of 30"-39" 
and wavelength ranges of 1355.08A-1355.43A and 
2794.10A-2794.45A, respectively. The time his¬ 
tory of the Mg II 2798. 8A intensity (integrated 
over 2798.65A-2799.00A, NUV continuum sub¬ 
tracted) was also calculated over the same slit 
positions. The continua and Mg II light curves 
appear to be similar. Figure [l] is a snapshot of this 
movie. Movie m2.mp4 shows the time evolution of 
the Fe xxi 1354.08A and Si IV 1402.77A spectra 
at the full-cadence (^9.5 s) in a smaller section 
of the slit. Note that different images have been 
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Fig. 1.— (A): IRIS/SJI 1330A image taken at 16:55:35 UT on 2014 September 6. (B)-(D): IRIS detector 
images taken at the same time in three spectral windows. (E) SDO/AIA 131A image taken at 16:55:32 
UT. (F): Light curves (arbitrarily scaled) of the GOES 1-8A flux and its time derivative, continua and 
Mg II 2798.8A intensities averaged over the slit positions of SO^-SO", and RHESSI 25-50 keV flux. In panels 
(A) and (E), the slit location is indicated by the white vertical line, and contours of RHESSI 25-50 keV 
X-Ray flux are overplotted on the images. The vertical line in panel (F) indicates the observation time of 
the images. Two online movies ml.mp4 and m2.mp4 are associated with this figure. 
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scaled differently in this movie to reveal a bet¬ 
ter contrast between the Fe xxi emission and the 
other emission features. 


2.2. Results and discussion 


Figure [T] shows that the IRIS slit crosses the 
edge of a ribbon, where one footpoint of the 
flare loops is anchored. The higher part of the 
loops is also covered by the slit. From the on¬ 
line movie ml.mp4 we can see that discernible 
Fe xxi 1354.08A emission first appears at about 
16:51:57 UT around the slit position of 28". 
About one minute later, stronger and entirely 
blueshifted Fe xxi emission appears around slit 
position 30", where a localized brightening can 
be identified in both the 1330A and 13lA im¬ 
ages. The entirely blueshifted Fe xxi line with¬ 
out a stationary component ha s been recent ly re¬ 
ported in seve r al other flares (Young et al. _ 20151: 
Li et al.l 20151: Polito et al. 20151 : ISadvkov et al 


2015 


Tian et all 2014al ). and is in agreement with 


hydrodynamic simulations of single flare loop (e.g., 
Emslie fe Alexander 1987 : Lin et al. 2009ajh En¬ 
hancement of both the FUV and NUV continua 
can be seen at slit positions 30"-39" after 16:53:51 
UT, ^2 minutes after the first Fe xxi 1354.08A sig¬ 
nature. During the impulsive phase of the flare, 
intensities of the 1330A and 13lA passbands, 
the transition region and chromospheric emission 
lines, and the FUV and NUV continua all increase 
dramatically. These intensities reach their max¬ 
ima at around 16:55:25 UT, almost the same as 
the peak times of the RHESSI HXR and GOES 
derivative. The 131A emission is likely domi- 
nated by Fe vui lines early in the impulsive phase 
( Boerner et al. 11 2012L iBrosius V Holmanll2012f ). In 


the impulsive phase of the flare, the Fe xxi emis¬ 
sion is found mainly at or close to the ribbon where 
the continua are enhanced, and the line is clearly 
blue shifted. A visual inspection of the movie also 
suggests that the Fe xxi blue shift reaches the 
maximum around the peak times of the RHESSI 
HXR and GOES derivative. The flare loops have 
little Fe xxi 1354.08A emission in the early stage 
of the impulsive phase. However, more and more 
Fe xxi emission appears to move upward from the 
loop footpoint to higher part in the later stage of 
the impulsive phase, and eventually fills the whole 
flare loops around the peak time of the GOES 
flux. This process provides strong support to the 


argument that the SXR emission in flare loops 
is produced through chromospheric evaporation, 
although we can not rule out the possibility that 
some of t he cor ona l plas ma might be heated in 
situ (e.g., iLiu et all 1201311 ). The Fe xxi emission 
from these flare loops is very strong and shows 
almost no Doppler shift. 

Figure [2]presents typical spectra at the flare rib¬ 
bons and flare loops in this flare. The ribbon spec¬ 
tra and loop spectra appear to be largely different. 
First, almost all emission lines as well as the FUV 
and NUV continua are much stronger at the rib¬ 
bons than in the loops. The Fe xxi 1354.08A line, 
however, is often stronger in the loops. This is be¬ 
cause hot plasma accumulates in the flare loops as 
the evaporation proceeds. Second, many narrow 
lines, which are normally invisible, now appear 
in the ribbon spectra. Most of these lines are 
singly ionized or neutral lines, and some of them 
have been identified with previous instruments 
( Cohenl 198ll: Sandlin et all 1986; Feld man et al 


Ekberg fe Feldman 2003 : Curdt et al.l 200 ll 
Tian et al.ll2009h.~Some l i nes an d their vac- 


199 


200 


mim wavelengths ( Kelly 1987 . 1979h have been 


marked in Figures [2] Third, the NUV spectra at 
the flare ribbons are domin ated by emission lines 
(see also lCheng et al.l2015b ). This is the case even 
at longer wavelengths around 2832A (not shown 
here), indicating photospheric heating through 
possi bly radiative backwarming (e.g., IXn et al. 1 


120041 ). The Mg II 2798.809A line, which is actu¬ 
ally a blend of two Mg II lines at 2798.754A and 
2798.822A, is the strongest emission line in the 
NUV spectrum except for the Mg II k and h lines. 
These two lines (Mg II 2798.8A line thereafter) 
are of particular interest as they can serve as a 
good diagnost ic of heating in the lower chromo¬ 
sphere ([Pereira et al.ll2015h . Away from the flare 
ribbons, the NUV spectra are similar to quiet Sun 
spectra with all lines except the Mg II k and h in 
absorption. Fourth, transition region lines such 
as C II and Si IV are often obviously broader and 
redshifted at the ribbons. 


In the following we focus on the flare rib¬ 
bon region (loop footpoint) crossed by the IRIS 
slit. Figure [3] shows the wavelength-time plots for 
Fe xxi 1354.08A and Si IV 1402.77A at four differ¬ 
ent positions of the ribbon. At slit positions 30.7" 
and 31.4", Fe xxi emission with a blue shift of 
^100 km s -1 are clearly present before the FUV 
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Fig. 2.— Typical IRIS spectra in the Ml.l flare on 2014 September 6. (A) Detector images taken at 16:56:13 UT. (B) Detector images 
taken at 17:08:59 UT. (C) Typical spectral line profiles at the flare ribbon (averaged over 10 pixels around slit position 34.9" in panel (A)) 
and flare loop (averaged over 10 pixels around slit position 16.6" in panel (B)). The flare ribbon spectra in the two Si IV spectral windows 
have been divided by 3 for the purpose of illustration. 
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Fig. 3.— Time evolution of the spectral line profiles of Fe xxi 1354.08A and Si IV 1402.77A at four 
different positions on the slit in the 2014 September 6 observation. The zero velocity has been set to the 
rest wavelength of the Fe xxi (upper panels) or Si IV (lower panels) line. Negative and positive velocities 
represent wavelengths shorter (blue shift) and longer (red shift) than the rest wavelength, respectively. The 
RHESSI 25-50 keV flux (purple), GOES 1-8A flux (white solid) and its time derivative (white dashed) are 
overplotted. Light curves of the UV continua (blue for FUV and green for NUV) and Mg II 2798.8A (cyan) 
intensities at the corresponding slit positions are also overplotted. 






continuum enhancement. This Fe xxi emission is 
probably from a loop with a footpoint not under 
the slit. At slit positions 33.4" and 34.7", there is 
essentially no Fe xxi emission prior to the contin¬ 
uum enhancement. The maximum continuum en¬ 
hancement occurs at or slightly after the peak time 
of the RHESSI HXR and GOES derivative (for the 
different spatial scales from which the X-Ray and 
UV curves are derived, we refer to a discussion 
later in this section). Since the UV continuum is a 
marker of chromospheric heating, this timing sug¬ 
gests that nonthermal electrons deposit energy in 
the chromosphere and heat the plasma there. Fol¬ 
lowing the maximum continuum enhancement, we 
see the appearance of largely blueshifted Fe xxi, 
especially at slit position 33.4". The blue shift of 
Fe XXI smoothly decreases to nearly zero within 
~3 minutes, whereas the HXR and UV continuum 
appear to decay even faster. 

Although a few profiles of Si IV 1402.77A are 
saturated around the times of maximum contin¬ 
uum enhancement, it is still clear that the inten¬ 
sities of Si IV, UV continua and Mg II 2798.8A 
all enhance at the same time and evolve in a 
very similar way, which suggests that the Si IV 
emission is responsible for the enhanced emission 
in lower layers and thus supports the idea that 
the FUV continuum enhancement at the flare rib¬ 
bons is driven by back-warming (heating from the 
downward emission component) from the strong 
C II and Si IV lines ( Dovle fe Phillips!Il992 ). An¬ 
other interesting feature of the Si IV line pro¬ 
files is the obvious asymmetry of the line profiles. 
Red wing asymmetries have been previously re¬ 
ported in the H<a line profiles acquired through 
some gro und based observations of flares (e.g., 
Ding et al. 1995 ). whereas observations by CDS 


and EIS usually reveal a nearly symmetric pro¬ 
files of the chromospheric and transition region 
lines during flares. The clear enhancement at the 
red wing suggests the presence of downflows with 
speeds of generally a few tens of km s -1 . Some of 
these downflows, e.g., those at the time of max¬ 
imum continuum enhancement, most likely indi¬ 
cate the ” chromo spheric condensation ” in explo¬ 


sive evaporation ( Fisher et al. 1985al) . The red 


wing enhancement can extend to more than 100 
km s -1 , suggestive of downflows at speeds greater 
than previously thought. Note that in many pre¬ 
vious observations the downflow speeds were es¬ 


timated based on a simple Gaussian fit to the 
line profiles of cool lines. The single Gaussian fit 
tends to underestimate the magnitude of red shift 
if strong red wing asymmetry is present. The red 
wing enhancement experiences a fast decay before 
reaching a stable level. The initial decay appears 
to be faster than the Fe xxi blue shift, consis¬ 
tent with the model of iFisher et al.l ( 19891 ) which 
predicts a chromospheric condensation time of ~1 
minute. We find that the red wing asymmetry 
is present even after the impulsive phase. These 
asymmetries indicate downflows with speeds of 
~30 km s -1 , which may be the ’’warm rain” result¬ 
in g from cooling o f the heated plasma, as discussed 
by 


Brosins 


We derive the Fe xxi 1354.08A line param¬ 
eters (peak intensity, Doppler shift and line 
width) in the evaporation flow through multi- 
Gaussian fit to the Fe xxi 1354.08A line and 
several lines nearby. In the spectral range 
of 1352.6A-1354.5A, the other major emission 
lines besides Fe XXI 1354.08A are the following: 
Si ii 1352.635A and 1353.718A, C 1 1352.745A and 
1354.288A, Fe n 1353.023A and 1354.013A. These 


six lines are all emitted from neutrals or singly 
ionized ions, thus are usually much narrower com¬ 
pared to the Fe xxi 1354.08A line. These lines also 
have zero or small Doppler shift (often red shift). 
The most prominent and blueshifted Fe xxi emis¬ 
sion is found around slit position 33.4" after the 
HXR peak. Before the HXR peak the Fe xxi emis¬ 
sion is found slightly to the south of this position 
(first around slit position 30.7" and then 31.4", 
see Figure [3]), which might be related to some pre¬ 
flare activity or some energy release early in the 
flare. We thus choose three slightly different slit 
positions during three time periods: slit position 
30.7" during 16:52:45-16:54:38 UT, 31.4" during 
16:54:48-16:55:07 UT, and 33.4" after 16:55:16 
UT. For each exposure, we average the line pro¬ 
files over 5 pixels around the corresponding slit 
position to improve the signal to noise ratio. A 
seven-component Gaussian function has been fit¬ 
ted to each of these spatially averaged line profiles. 
Figure |H shows three examples of multi-Gaussian 
fit. In some spectra there appear to be some 
additional emission peaks (e.g., Figure ID^C)), sug¬ 
gesting the presence of some unidentified lines. 
Most of these lines are also very narrow and often 
weak in intensity, and their presence should not 
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Wavelength (Angstrom) 

Fig. 4. — Three examples of multi-Gaussian fit to Fe xxi 1354.08A and several lines nearby in the 2014 
September 6 observation. The observed and fitted line profiles are shown as the histograms and smooth 
curves, respectively. The dashed line indicates the rest wavelength of Fe xxi 1354.08A. 
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Wavelength (Angstrom) 

Fig. 5.— Line profiles of the 1354A (A and B) and 1349A (C) spectral windows at 16:55:16 UT and 16:55:25 
UT on 2014 September 6, at slit position 33.4". As a comparison the line profile at 16:57:57 UT is re-scaled 
and overplotted in panels (A) and (B). The dashed line indicates the rest wavelength of Fe xxi 1354.08A. 
Several previously unidentified lines are marked as the V followed by the approximate wavelengths. 
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significantly affect the derived line parameters of 
Fe xxi 1354.08A. 


At first sight the Fe xxi 1354.08A line seems 
absent at slit position 33.4" from 16:55:16 to 
16:55:25 UT. However, a careful examination of 
the spectra at 16:55:16 and 16:55:25 UT (Fig¬ 
ure 02) suggests some extra emission around the 
short wavelength edge of the spectral window. To 
demonstrate this we rescaled the line profile ob¬ 
tained at 16:57:57 UT so that the line free part 
of the spectrum from 1355.lA-1355.4A matches 
those at 16:55:16 and 16:55:25 UT. We see obvi¬ 
ous enhancement of the continuum level between 
the emission lines around the short wavelength 
edge of the spectral window at 16:55:16 UT and 
16:55:25 UT. Moreover, the C I 1352.7451 line is 
obviously weaker than the Fe n 1353.0231 line in 
all the selected spectra except the two obtained 
at 16:55:16 and 16:55:25 UT. It is likely that this 
bulk enhancement is due to the highly blueshifted 
(^200-300 km s -1 ) Fe xxi line. However, we 
realize that it may also be complicated by the 
redshifted emission of the cool lines. For instance, 
the continuum enhancement around 1353.2A may 
be contributed by the enhanced red wing em is¬ 


sion o f the Fe II 1353.023A line. lYoung et al 


(2015) noticed that several Si II lines sometimes 
show broad and extended red wings. This is 
also the case in our data. Strong and broad 
redshifted emission features are clearly present 
in the Si II 1353.718A line profiles at 16:55:16 
and 16:55:25 UT (Figure [5f A)-(B)). Similar fea¬ 
tures are also found for another two Si II lines at 
1348.543A and 1350.057A (Figure [5fC)), demon¬ 
strating that they are not unidentified lines. Thus, 
the enhanced emission of C I 1352.745 A and 
the continuum around 1352. 9A may actually be 
caused by the broad redshifted component of the 
Si II 1352.635A line. Because of this complication, 
we conclude that highly blueshifted Fe XXI line 
possibly exists but cannot be unambiguously iden¬ 
tified at 16:55:16 and 16:55:25 UT. 


The derived Fe xxi 1354. 08 A line parameters 
are plotted in Figure [6l Our study is among one 
of the few studies where hot evaporation flows 
are tracked at a cadence of ~10 s. Using CDS 
data, the blue shift of Fe xix 592. 23 A was t racked 
with a cadence of 9.8 s in tw o flares ( Brosiusl 
20031 : iBrosius fe Phillip i 12004 ). A cadence of 
7 s was used in the observation presented by 


Brosins fe Holmanl (|2010h but no blueshifts were 
found for the hot Fe xix line. With EIS observa¬ 


tions, iBrosiusI ( 201 3) reported the change of the 


intensity and Doppler shift of Fe xxm 263.8A in 
a Cl flare with a cadence of 11s. We find that the 
line width values (1/e width) are mostly in the 
range of 60-90 km s -1 . Given a thermal broaden¬ 
ing of 57.8 km s -1 (under ionization equilibrium) 
and an instrumental broadening of ^4 km s -1 
(see the discussion i n Section S5 of th e Supple¬ 
mentary Materials of lTian et al.l[2Q14b )« the non- 
thermal broadening has been found to be in the 
range of 15-70 km s -1 . The line intensity ap¬ 
pears to peak around the HXR peak time. We 
have to bear in mind that the line parameters 
are measured at the loop footpoint and that the 
Fe xxi emitting plasma continuously moves to 
higher part of the flare loops. As we explained 
above, the maximum blue shift of Fe xxi ap¬ 
pears in the time range of 16:55:16 UT - 16:55:35 
UT, roughly the same as the peak times of the 
RHESSI HXR and the derivative of GOES SXR. 
The blue shift quickly decreases to ^20 km s -1 
within ^3 minutes, before the start (17:02 UT) 
of the second episode of energy release in this 
flare. The uncertainty on the rest wavelength 
of the Fe xxi line may slightly change the de¬ 
rived Doppler shifts. In our calcu latio n we use 
1354.08A as the rest wavelength (|Sandlin et al.l 


1977). Usi ng a wavel ength of 1354.064A given by 
Feldman et al. feOOOl) will lower the derived blue 


shifts by ^3.5 km/s. While using a wav elength 
of 1354.106A found by Young et al. ( 20151 ) will in¬ 
crease the derived blue shifts by ^5.8 km/s. The 
continua and Mg II 2798.8A intensities averaged 
over the whole covered ribbon (slit positions of 
30"-39") reach their maxima at around 16:55:25 
UT, which is almost the same time of the max¬ 
imum blue shift. Note that the Fe xxi emission 
comes from relatively higher part of the loop legs, 
whilst the continua and Mg II emission come from 
lower part of the loop legs. Ideally one should 
compare the Doppler shift of Fe xxi and the con- 
tinua/Mg II intensities in the same loop. However, 
it is difficult to do do due to the line-of-sight su¬ 
perposition effect. 


We have to mention that although the RHESSI 
HXR and GOES SXR are both integrated over 
the whole Sun, a comparison of them with the UV 
emission and Fe xxi 1354.08A line parameters at 
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Start Time (06-Sep-14 16:52:30) 


Fig. 6.— Time history of the line parameters of Fe xxi 1354.08A, the GOES 1-8A flux and its time derivative, 
RHESSI 25-50 keV flux, continua and Mg II 2798.8A intensities averaged over slit positions of SO^-SQ", in 
the 2014 September 6 observation. The error bars represent the l-cr uncertainties from the multi-Gaussian 
fitting. Note that the Fe xxi line cannot be unambiguously identified at 16:55:16 UT and 16:55:25 UT, so 
there are no measurements of the Fe xxi line parameters at these two times. The black dashed line in panel 
(B) is the exponential fit of the time evolution of Fe xxi blue shift. The variable t is the time (second ) after 
16:55:35 UT. 
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the ribbon is still meaningful. First, most of the 
full-Sun RHESSI HXR flux is emitted from the 
two conjugate footpoints (see Figure [T}, each lo¬ 
cated on one of the two ribbons. Connected by the 
same flaring loop, such HXR footpoints have lig ht 


curves of the same temporal evolution (Liu et al 


2009bf ). Therefore, the full-Sun HXR flux 


is pro¬ 
portional to that of the footpoint crossed by the 
IRIS slit. The Neupert effect suggests that the 
GOES derivative is a proxy for the variation of 
the HXR flux, so it can be compared with the 
Fe xxi line parameters at the ribbon as well. Sec¬ 
ond, a recent study clearly reveals that the tem¬ 
poral evolution of evaporation flows at different 
spa tial pixels wit hin the ribbon is strikingly simi¬ 
lar (iGraham fc Gauzzl2015l ), although we keep in 
mind that the maximum flow speed might occur 
at different times at different pixels in some flares. 


To quantify possible correlation of the Fe xxi blue 
shift with other parameters, we present various 
scatter plots in Figure 0 for measurements at slit 
position 33.4" after 16:55:30 UT. Larger blue shifts 
are found when the intensities of Mg II 2798.8A 
and FUV/NUV continuum are larger, supporting 
the scenario of evaporation flows resulting from 
heating in the lower chromosphere. 


The Fe xxi blue shift also correlates with the 
RHESSI HXR and GOES derivative, meaning that 
the evaporation flow speed depends on the rate of 
energy input at the flare ribbon. Perhaps the most 
intriguing correlation is the one between the blue 
shift and GOES derivative, which is nearly lin¬ 
ear with a correlation coefficient of 0.98. In the 
thick-target flare model, most of the accelerated 
electrons produce HXR through bremsstrahlung 
in the chromosphere. The HXR flux is positively 
correlated with the instantaneous energy deposi¬ 
tion rate by the nonthermal electrons precipitat¬ 
ing to the chromosphere. The electron energy de¬ 
position may be lost through radiation, or go to 
the heating of the chromospheric plasma. Our 
results suggest that the evaporation flow is en¬ 
hanced when the energy deposition rate is higher, 
thus indicating that a significant fraction of the 
electron energy goes into chromospheric heating 
and the resultant evaporation in this flare. The 
evaporation flow velocity a ppears to dep e nd on 
the heating rate. Recently iKrucker et al ] (I 2015 I) 
claimed that energy in high-energy electrons (>30 
keV) is not responsible for chromospheric evapo¬ 


ration of hot plasma but is lost through radiation 
in th e optical range. T hrough analysis of an XI 
flare, iHeinzel fe Kleintl ( 20141 ) concluded that the 
enhanced NUV continuum is the hydrogen recom¬ 
bination Balmer-continuum emission following an 
impulsive ionization caused by the precipitation of 
electrons in the chromosphere, which also suggests 
the loss of electron energy through radiation. Our 
results are not necessarily in contradiction with 
these previously reported results, as the amount of 
chromospheric heating and optical radiation may 
both increase as the deposited energy increases. 

It is interesting that the Fe xxi blue shift is bet¬ 
ter correlated with the SXR derivative than the 
HXR and UV continua, although there may be 
some uncertainty due to the different spatial scales 
of integration. We offer possible explanations for 
this behavior: The blue shift represents how fast 
the heated plasma is transported into the flare 
loops, which may be more directly related to the 
energy deposition rate into the flare loops, as rep¬ 
resented by the SXR derivative. The HXR and UV 
continua are also correlated with the energy depo¬ 
sition rate, but they are, however, produced at the 
loop foopoints/ribbons. Therefore a positive cor¬ 
relation between them and the Fe xxi blue shift is 
expected, but may not be very tight, especially 
when we consider the nonlinearity in the com¬ 
plex processes involving Coulomb collisions, ion¬ 
ization, free-bound and bremsstrahlung radiation, 
and hydrodynamic response of the atmosphere. 
Another possible explanation is that the chro¬ 
mospheric evaporation is mainly driven by ther¬ 
mal c onduction rather than direct electron heatin g 
(e.g., IZarro fe Lemenl[l988[ iBattaedia et al.l[2QQ9h . 


which may lead to a better correlation between the 
blue shift and SXR derivative. 

Although there is a weak trend that smaller 
blue shifts are associated with smaller nonthermal 
widths, we do not see a strong correlation between 
th e blue shift and no nthermal width as reported 
bv IPolito et al. (|2Q15| ) in another flare. 


3. The XI. 6 flare on 2014 September 10 
3.1. Observations and data reduction 

The second flare we have analyzed is the XI.6 
flare peaking at about 17:45 UT on 2014 Septem¬ 
ber 10 from the sigmoidal region in NO A A AR 
12158. The IRIS observation lasted from 11:28 
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Fig. 7.— Scatter plots showing the correlation between Fe xxi blue shift and other parameters in the Ml.l 
flare on 2014 September 6. The linear Pearson correlation coefficient (C.C.) is marked in each panel. Only 
data points measured after 16:55:30 UT are used for the scatter plots. 
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Fig. 8.— Context images of the XI.6 flare on 2014 September 10. (A)-(B) Still images of AIA 1600A and 
131A from an online movie m3.mp4 showing the evolution of the entire flare region. (C)-(H): IRIS/S JI 1400A, 
AIA 131A and XRT ALthick images of the eastern ribbon taken around 17:29:47 UT and 17:31:56 UT. The 
red arrow in each panel indicates a fast apparent motion along the ribbon. The vertical line indicates the 
slit location. The box shown in panels (E) and (H) marks the region where the XRT intensity is averaged. 
(I)-(K): Space-time maps for the cut indicated by the red arrow in panels (C)-(H). The slope of the black 
line is used to estimate the speed of the apparent motion. An online movie m4.mp4 shows the evolution of 
the XRT intensity in which the FOV of Figure |H(C)-(H) is indicated by the rectangle. 
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Fig. 9.— (A): IRIS/SJI 1400A image taken at 17:34:11 UT on 2014 September 10. (B)-(D): IRIS detector 
images taken at the same time in three spectral windows. (E) SDO/AIA 131A image taken at 17:34:11 
UT. (F): Light curves (arbitrarily scaled) of the GOES 1-8A flux and its time derivative, continua and 
Mg II 2798.8A intensities averaged over 5 pixels around the slit position 117.8", XRT ALthick intensity 
integrated over the box shown in Figure [HfE) and its time derivative. In panels (A) and (E), the slit location 
is indicated by the white vertical line. The vertical line in panel (F) indicates the observation time of the 
images. Two online movies m5.mp4 and m6.mp4 are associated with this figure. 
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UT to 17:59 UT, which covers the whole impul¬ 
sive phase and part of the decay phase of the flare. 
The IRIS slit was oriented in the nominal north- 
south direction. With a pointing of (-168", 125"), 
the slit crossed two locations of the eastern rib¬ 
bon in this flare (Figure [8]). The observing mode 
was also sit-and-stare. The spatial pixel size, FUV 
spectral pixel size, NUV spectral pixel size, and 
cadence were 0.166", ~0.026 A, ~0.025 A, and 
~9.5 seconds, respectively. Slit-jaw images (SJI) 
were taken at a cadence of ^19 seconds in each 
of the 1400A and 2796A passbands. The start¬ 
ing exposure time was ^8 seconds, and the au¬ 
tomatic exposure control successfully reduced the 
exposure time of the NUV spectra to 2.4 seconds 
after 17:27:09 UT. during the flare. The level 2 
data are used here. The absolute wavelength cal¬ 
ibration and cross-channel coalignment have been 
done similarly to the 2014 September 6 dataset. 
We notice that the Si IV spectra obtained by IRIS 
at the southern part of the eastern ribbon have 
been analyzed bv lLi fc Zhand ( 2015 ) for the study 
of slipping magnetic reconnection. The IRIS data 
i n this event has also been used by ICheng et al.l 
( 2015aj ) and iLi et al.l ( 2015f ) to study the forma¬ 
tion of flux ropes and quasi-periodic pulsations, 
respectively. In this paper we focus mainly on the 
Fe xxi 1354.08A line emission at the northern part 
of the eastern ribbon. 


The AIA images have also been analyzed. Sim¬ 
ilar to the 2014 September 6 dataset, we used 
the AIA 1600A images for the coalignment be¬ 
tween IRIS and AIA data. For the AIA data we 
mainly focus on the 13lA image sequence, which 
has a cadence of 12 seconds and a pixel size of 
0.6". We notice that many pixels are saturated 
in the 13lA images taken around the flare peak 
time. An online movie m3.mp4 shows the evolu¬ 
tion of the entire flare region in the AIA 1600A and 
131A passbands. 

RHESSI data are not available for this flare. 
SXR images were taken with the HINODE/XRT 
instrument in several filters, although the field 
of view (FOV) of XRT only has partial overlap 
with that of IRIS/SJI. XRT images taken in the 
ALthick, Be_thin and Be_med filters appear to be 
similar. We thus use only the ALthick images 
in this paper. To coalign these images we have 
performed a cross-correlation of every 5 images. 
These internally aligned XRT images have then 


been aligned with the AIA 131A images, and thus 
also the IRIS images. A movie of the XRT images 
is available online (m4.mp4). 

Figure [8j A)-(F) presents the IRIS/SJI 1400A, 
AIA 131A and XRT ALthick images taken at two 
different times. A dashed line was drawn to in¬ 
dicate a fast motion seen in these images. The 
space-time maps for this dashed line are shown in 
Figure |H(G)-(I). 

We present the time sequence of IRIS obser¬ 
vation (including SJI 1400A images and three 
spectral windows) and the corresponding AIA 
13lA images in an online movie m5.mp4. The 
GOES 1-8A flux and XRT ALthick flux integrated 
over the box marked in Figure [8fE)&(H) are also 
shown in the movie. These two light curves appear 
to have a similar trend, and the general trend of 
their time derivatives is also similar. We have no¬ 
ticed an interesting Fe xxi 1354.08A emission fea¬ 
ture around the slit position 117.8". Since our fol¬ 
lowing analysis is focused mainly on this feature, 
we calculated the NUV and FUV continua by aver¬ 
aging the IRIS spectra over 5 pixels around slit po¬ 
sition 117.8" and wavelength ranges of 1355.08A- 
1355.43A and 2794.10A-2794.45A, respectively. 
The time history of the Mg II 2798.8A intensity 
(integrated over 2798.65A-2799.00A, NUV con¬ 
tinuum subtracted) was also calculated over the 
same slit positions. Figure [9] is a snapshot of this 
movie. Movie m6.mp4 shows the time evolution of 
the Fe xxi 1354.08A and Si IV 1402.77A spectra 
at the full-cadence (~9.5 s) in a smaller section 
of the slit. Similar to movie m2.mp4, different 
images have been scaled differently in m6.mp4 to 
reveal a better contrast between the Fe xxi emis¬ 
sion and the other emission features. 

3.2. Results and discussion 

From the XRT ALthick and AIA 131A images 
we can see that the sigmoid consists of two “J”- 
shape loops which evolve rapidly during the flare. 
The IRIS slit crosses the ribbon where the east¬ 
ern “J” loop is rooted. A fast motion of bright 
front can be clearly identified from 17:29:48 UT 
to 17:32:37 UT at the northern part of this ribbon 
(Figure [8}. This apparent motion is seen in not 
only the hot XRT ALthick and AIA 131A chan¬ 
nels, but also in the IRIS/SJI 1400A passband 
which samples mainly the Si IV lines formed at 
a temperature of ^10 4-9 K. The similar dynamics 
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of both the cool and hot plasma is consistent with 
the scenario of chromospheric heating and evapo¬ 
ration in sequentially reconnected loops, thus sup¬ 
po rting the interpreta tion of slipping reconnection 
by III fe Zhang ( 201 5). The apparent speed of the 
motion has been estimated to be ^200 km s -1 
from the space-time maps (Figure[E(G)-(I)), which 
is much higher than the sound speed at a temper¬ 
ature of ^10 4 ' 9 K. We believe that the fast ap¬ 
parent motion is not real mass motion but just a 
manifestation of the progression of the magnetic 
reconnection site along the flare arcade. In other 
words, the fast apparent motion is ultimately due 
to a succession of chromospheric footpoints being 
heated with nonthermal electrons accelerated by 
reconnection in a succession of loops as the mag¬ 
netic reconnection site progresses along the flare 
arcade. 

Figure [9] shows clear enhancement of both the 
FUV and NUV continua at both the southern and 
northern parts of the covered ribbon. In the south¬ 
ern part, the Mg II 2798.8A line turns from absorp¬ 
tion to emission at about 17:24:29 UT and clear 
continuun enhancement can be identified ^20 s 
later (Movie m5.mp4). Obvious red wing enhance¬ 
ment frequently appears in the line profiles of C I, 
Si IV and Mg II, indicating downward motion of 
the chromosphere at a few tens km s -1 due to the 
heating related overpressure. However, no obvi¬ 
ous signature of Fe xxi 1354.08A is found at this 
part, suggesting that the chromospheric heating 
is likely not intense enough to produce ~11 MK 
plasma there. 

The chromospheric heating appears much more 
efficient in the northern part of the ribbon, where 
we see greatly enhanced continua and largely 
blueshifted Fe xxi 1354.08A. From the online 
movie m5.mp4 we can see that the continuum 
enhancement and Mg II 2798.8A emission first ap¬ 
pears around 17:25:25 UT at slit position ~126". 
As the ribbon extends in the north-south direction 
after 17:30:07 UT, the slit positions where the con¬ 
tinua and Mg II 2798.8A intensities are enhanced 
expand and cover the range of 110"-134" around 
the flare peak time. Clear Fe xxi 1354.08A sig¬ 
nature can be identified at slit position ~128" at 
17:28:52 UT. The blueshifted Fe xxi emission can 
be more easily identified as the ribbon extends. 
The largest blue shift appears to be found at the 
southern edge of the ribbon, where the continua 


and Mg II 2798. 8A intensities are also the largest. 
This suggests that the largest blue shift is likely 
associated with the newly heated plasma. Be¬ 
fore the GOES SXR peak, we can clearly see that 
entirely blueshifted Fe xxi emission features con¬ 
tinuously migrate to the rest wavelength position 
of the line, leading to more and more nearly sta¬ 
tionary Fe xxi emission at the locations of the 
hot flare loops. This behavior provides support 
to the scenario of the hot emission in flare loops 
being produced through an accumulation of the 
evaporated plasma. A smooth transition from the 
largest blue shift at the outer edge of the ribbon 
to a zero shift in the f l are lo ops, which has been 
found by I Young et al. (|2015h in another flare, is 
clearly present in the later stage of the impulsive 
phase (^17:35 UT - 17:45 UT). The main evap¬ 
oration process appears to cease after the GOES 
SXR peak, when the blue shift of Fe xxi almost 
disappears and the intensity of the nearly station¬ 
ary Fe xxi emission decreases with time due to 
the cooling of the hot plasma. 

Figure [TO] shows the wavelength-time plots for 
Fe xxi 1354.08A and Si iv 1402.77A at four differ¬ 
ent positions in the northern part of the ribbon. 
Largely blueshifted Fe xxi emission appears at or 
after the time of sudden enhancement of the FUV 
continuum at all positions, again supporting the 
scenario of evaporation flows driven by heating in 
the lower chromosphere. The maximum blue shift 
of Fe xxi is on the order of ~250 km s -1 at the 
slit positions of 114.3", 117.8" and 120.6". The 
blue shift smoothly decreases to nearly zero within 
~9 min, before the peak of the GOES flux, at all 
four positions. Note that the appearance of con¬ 
tinuum enhancement and Fe xxi emission at dif¬ 
ferent times at the four different positions is just 
a reflection of the ribbon expansion to the south. 

Although saturated around the times of max¬ 
imum continuum enhancement, the Si IV spectra 
presented in Figure [TOl EWH) reveal some results 
similar to those in the 2010 September 6 event. 
The strong correlation of Si IV 1402.77A, continua 
and Mg II 2798.8 A line intensities supports the 
suggestion of FUV continuum enhancement being 
driven by back-warming from the C II and Si IV 
lines ( Dovle fe Phillipsl Il992h . The obvious red 
wing asymmetry in the impulsive phase and decay 
phase are like ly signatures of the chromospheric 
condensation ( Fisher et al.l Il985ah and cooling of 
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Fig. 10.— Time evolution of the spectral line profiles of Fe xxi 1354.08A and Si IV 1402 .Ilk at four 
different positions on the slit in the 2014 September 10 observation. The GOES 1-8A flux (white solid) and 
its time derivative (white dashed), XRT ALthick intensity integrated over the box shown in Figure [8fE) 
(purple solid) and its time derivative (purple dashed), continua (blue for FUV and green for NUV) and 
Mg II 2798.8A (cyan) intensities at the corresponding slit positions are overplotted. 
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Fig. 11.— Three examples of multi-Gaussian fit to Fe xxi 1354.08A and several lines nearby in the 2014 
September 10 observation. The observed and fitted line profiles are shown as the histograms and smooth 
curves, respectively. The dashed line indicates the rest wavelength of Fe xxi 1354.08A. 


21 



































Wavelength (Angstrom) 

Fig. 12. — Line profiles of the 1354A (A and B) and 1349A (C) spectral windows at 17:30:45 UT and 
17:31:41 UT on 2014 September 10, at slit position 117.8". As a comparison the line profile at 17:41:33 
UT is re-scaled and overplotted in panels (A) and (B). The dashed line indicates the rest wavelength of 
Fe xxi 1354.08A. Several previously unidentified lines are marked as the V followed by the approximate 
wavelengths. 
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the heated plasma ( Brosiusl 2003| ? respectively. 


Since the ribbon moves a lot in the slit direc¬ 
tion and there are many loops overlapping with 
each other, it is difficult to identify and track 
the Fe xxi emission in a single loop. Fortu¬ 
nately, we find that the Fe xxi emission fea¬ 
ture at the slit position 117.8" appears not mov¬ 
ing in the slit direction from ^17:32:28 UT to 
17:41:23 UT. We track the time evolution of the 
Fe xxi line parameters at this position. Simi¬ 
lar to the 2010 September 6 event, we average 
the line profiles over 5 pixels around slit posi¬ 
tion 117.8" and derive the Fe xxi 1354.08A line 
parameters through a seven-component Gaussian 
fit to the Fe xxi 1354.08A line and six cool lines 
nearby. Figure fill shows three examples of multi- 
Gaussian fits. We notice that extra emission 
around the short wavelength edge of the spectral 
window, similar to those seen in the 2014 Septem¬ 
ber 6 flare (Figure [5j), is also clearly present in 
the line profiles obtained from 17:30:35 UT and 
17:32:19 UT. Figure fl2f A)-(B) show two exam¬ 
ples. There the line profile obtained at 17:41:33 
UT has been re-scaled to let the line free part 
of the spectrum from 1355.lA-1355.4A matches 
those at 17:30:45 UT and 17:31:41 UT. The extra 
emission may indicate a very large blue shift of 
Fe xxi, ~300 km s -1 or even greater since part of 
the Fe xxi line profile may have been shifted out 
of the spectral window. However, as we discussed 
above, the presence of a broad redshifted compo¬ 
nent (e.g., at 17:30:45 UT) or red wing enhance¬ 
ment (e.g-, at 17:31:41 UT) in three other Si II lines 
(1348.543A and 1350.057A shown in FigurefT2fC). 
1353.718A shown in Figure fl2f AU(B)) suggests a 
significant contribution of the extra emission by 
the redshifted Si II 1352.635 A line. We decide not 
to perform multi-Gaussian fitting to the line pro¬ 
files obtained from 17:30:35 UT and 17:32:19 UT 
since it is difficult to decompose the contribution 
from the highly blueshifted Fe xxi line and the 
redshifted Si II emission. 


The derived Fe xxi 1354.08A line parameters 
and the 1 -a uncertainties are plotted in Figure fl3l 
The 1 /e line width values are mostly in the range 
of 60-80 km s -1 , leading to a nonthermal broad¬ 
ening of 10-50 km s -1 after subtracting the ther¬ 
mal and instrumental broadening. Different from 
the 2014 September 6 event, the line intensity de¬ 
creases slightly before increasing continuously till 


the end of the measurements at 17:41:23 UT. As 
we explained above, Fe xxi emission with a blue 
shift larger than 200 km s -1 first appears in the 
time range of 17:30:35 - 17:32:28 UT, same as or 
slightly after the time of the sudden enhancement 
of the continua and Mg II 2798.8A intensities. 
The maximum blue shift appears to precede the 
peaks of the SXR (GOES and XRT) derivatives 
by ~0.5-2.0 minutes. This delay might be quali¬ 
tatively consistent with a recent radiative hydro- 
dynamic sim ulati on (RADYN) _of chromospheric 
evaporation ( Rubio da Costa et al ] 2015h . where 
the maximum upflow speed is achieved before the 
maximum of the injected electron flux, and thus 
HXR flux or time derivative of SXR flux. How¬ 
ever, we have to remember that the maximum blue 
shift actually occurs at different ti mes at differ¬ 
ent pixels within the ribbon ( Graham &; Gauzzil 
Note that the XRT flux is integrated 


20151 


over a small region around the position where the 
Fe xxi 1354.08A line parameters are measured. 


The blue shift smoothly decreases to ^0 km s -1 
in ~9 minutes. The time evolution of the blue 
shift from 17:32:28 to 17:41:23 UT appears to 
be well fitted with an exponential function. It 
suggests that the evaporation process evolves 
smoothly, perhaps indicating a continuous en¬ 
ergy input. However, there is an uncertainty 
as to whether overlying loops corresponding to 
different footpoints may be contributing to the 
Fe xxi 1354.08A emission feature at slit position 
117.8". If this was the case, then we would expect 
very broad Fe xxi line profiles or perhaps multiple 
components which does not appear to be the case. 


Figure [14] presents various scatter plots for the 
relationship between Fe xxi blue shift and other 
measured quantities. As expected from the sce¬ 
nario of evaporation flows being driven by chro¬ 
mospheric heating, larger blue shifts are found 
when the intensities of Mg II 2798.8A and the 
FUV and NUV continua are larger. The corre¬ 
lation appears to be linear for blue shifts larger 
than ^130 km s -1 . The chromospheric intensities 
do not drop too much when the Fe xxi blue shift 
reaches smaller values, simply because they drop 
faster than the decrease of the blue shift. Similar 
to the 2014 September 6 event, here we also find 
a good correlation between the Fe xxi blue shift 
and the SXR derivative. The best correlation is 
still the one between the Fe xxi blue shift and the 
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Fig. 13.— Time history of the line parameters of Fe xxi 1354.08A, the GOES 1-8A flux and its time 
derivative, XRT ALthick intensity integrated over the box shown in Figure [3(E) and its time derivative, 
continua and Mg II 2798.8A intensities averaged over 5 pixels around the slit position 117.8", in the 2014 
September 10 observation. The error bars represent the 1-cr uncertainties from the multi-Gaussian fitting. 
Note that highly blueshifted (^300 km s -1 or larger) Fe xxi line may be present from 17:30:35 and 17:32:19 
but cannot be unambiguously identified. The black dashed line in panel (B) is the exponential fit of the time 
evolution of Fe xxi blue shift. The variable t is the time (second ) after 17:32:28 UT. 
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Fig. 14.— Scatter plots showing the correlation between Fe xxi blue shift and other parameters in the XI.6 
flare on 2014 September 10. The linear Pearson correlation coefficient (C.C.) is marked in each panel. 
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derivative of GOES, which has a linear correlation 
coefficient of 0.94. As we explained above, this 
correlation means that the evaporation flow veloc¬ 
ity largely depends on the energy deposition rate, 
and thus the heating rate. 

The scatter plot for the relationship between 
the blue shift and nonthermal width (Figure fl4fF)) 
appears surprising. We see a positive correlation 
when the blue shift is larger than ^100 km s -1 
and a negative correlation when the blue shift is 
smaller than ~100 km s -1 . More case studies will 
tell how common this characteristic is. 


4. Summary 


Using IRIS observations, we have tracked the 
complete evolution of the hot evaporation flows 
continuously in the Ml.l flare occurring on 2014 
September 6 and the XI.6 flare on 2014 Septem¬ 
ber 10. For both flares, the hot evaporation flows 
marked by the blueshifted Fe xxi 1354.08A emis¬ 
sion are tracked with a cadence higher than 
10 seconds. The relationship between chromo¬ 
spheric evaporation and heating, as indicated 
by the enhancement of the UV continuum and 
Mg II 2798.8A line emission, has been studied as 
well. We have also investigated the correlation be¬ 
tween the evaporation flow velocity and the energy 
deposition rate, as represented by the HXR flux 
measured by RHESSI or the time derivative of the 
SXR flux observed by GOES and HINODE/XRT. 
The main results are summarized in the following. 

(1) The NUV spectra at the ribbons of both 
flares are dominated by emission lines. The is true 
for wavelengths at least between 279lA-2834A 
and indicates strong heating in both the chromo¬ 
sphere and photosphere. 

(2) The Fe xxi line is completely blueshifted at 
or near the flare ribbons indicated by the enhance¬ 
ment of the UV continuum. The absence of the 
rest component, which was often observed in flare 
lines by previous instruments, agrees with simula¬ 
tions of chromosphe ric evapo ration in a singln flare 


loop (e.g., lEmslie fe Alexanderl 119871 : iLiu et al. 


l2009aj ) and suggests that the resolution of IRIS 
is high enough to separate the evaporation flows 
from the surrounding hot plasma in flare loops. 

(3) The flare loops have little Fe xxi emis¬ 
sion in the early stage of the impulsive phase. 
Fe xxi emission appears to shift from the foot- 


points to the higher parts of the flare loops in the 
later stage of the impulsive phase, until the entire 
flare loop emits in Fe xxi at the time of the peak 
in the GOES flux. This behavior provides support 
to the scenario of the SXR emission in flare loops 
being produced through an accumulation of the 
evaporated hot plasma from the chromosphere, al¬ 
though we can not rule out the possibility of cer¬ 
tain degrees of in-situ heating in the corona. 

(4) The maximum blue shift of Fe xxi is on 
the order of 200 km s _1 , though possible signa¬ 
tures of ^300 km s _1 or larger blue shifts have also 
been identified. The maximum blue shift appears 
approximately at the same time as or slightly af¬ 
ter the impulsive enhancement of the FUV/NUV 
continuum and Mg II 2798.8A emission, indicating 
that the hot evaporation flow is closely related to 
intense heating of the lower chromosphere. 

(5) After reaching the maximum, the blue shift 
decreases smoothly to almost stationary within a 
few minutes. The smooth decrease may suggest 
that the energy input into the lower atmosphere 
is smooth and continuous in these flares. The ve¬ 
locity decrease appears to be exponential in time, 
especially for the XI.6 flare. 

(6) The Fe xxi blue shift appears to correlate 
with the hard X-Ray flux (or derivative of soft 
X-Ray flux), which is in general agreement with 
the picture of chromospheric evaporation driven 
by nonthermal electron heating. A nearly lin¬ 
ear correlation between the blue shift and GOES 
derivative is found in both flares, suggesting that 
the evaporation flow velocity largely depends on 
the rate of energy deposition or heating. 

(7) The intensities of the FUV continuum, 
NUV continuum, and the chromospheric and 
transition region lines appear to be correlated. 
This correlation supports the suggestion of FUV 
continuum enhancement in flares being driven 
by back-warming from the C II and Si IV lines 


(jDovle &; Phillips! ! 1992b The decay of these inten¬ 


sities is faster than the decrease of the Fe xxi blue 
shift. 


(8) The chromospheric and transition region 
lines often show obvious enhancement in the red 
wing of the line profiles at the flare ribbons. 
Several Si II lines sometimes even show a dis¬ 
tinct broad and redshifted component. The red 
wing enhancement and redshifted component sug- 


26 











gest that the speed of chromospheric condensa¬ 
tion might be larger (up to ^100 km s -1 ) than 
previously found from hydrodynamic simulations 
and from single Gaussian fits to the cool lines of 
EIS and CDS. Their presence also complicates the 
detection of highly blueshifted Fe xxi 1354.08A 
around the time of maximum continuum enhance¬ 
ment. The red wing enhancement is found in not 
only the impulsive phase but also the decay phase. 
The latter is likely caused by cooling of the heated 
plasma in the corona. 

(9) In the XI.6 flare, we see a bright intensity 
front propagating with an apparent speed of ^200 
km s -1 along the ribbon in IRIS/SJI 1400A, XRT 
ALthick and AIA 131A images. This fast apparent 
motion is not real mass motion but just a reflection 
of the chromospheric heating and evaporation in 
sequentially reconnected loops. 

IRIS is a NASA small explorer mission devel¬ 
oped and operated by LMSAL with mission oper¬ 
ations executed at NASA Ames Research center 
and major contributions to downlink communi¬ 
cations funded by the Norwegian Space Center 
(NSC, Norway) through an ESA PRODEX con¬ 
tract. Hinode is a Japanese mission developed 
and launched by ISAS/JAXA, with NAOJ as 
domestic partner and NASA and STFC (UK) 
as international partners. It is operated by 
these agencies in co-operation with ESA and 
the NSC (Norway). This work is supported by 
NASA contracts 8100002705, NNM07AB07C and 
SP02H1701R from LMSAL to SAO, NASA grants 
NNX13AE06G and NNX15AJ93G, and NSF grant 
AGS-1159353. H.T. thanks Fatima Rubio da 
Costa and Xin Cheng for insightful discussion. 
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